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We present a precise theoretical prediction for the decay width of the bound state of two electrons and a 
positron (a negative positronium ion), r(Ps^) = 2.087 085(12)/ns. We include O(q^) effects of hard virtual 
photons as well as soft corrections to the wave function and the decay amplitude. An outcome of a large-scale 
variational calculation, this is the first result for second-order corrections to a decay of a three-particle bound 
state. It will be tested experimentally in the new positronium-ion facility in Garching in Germany. 

PACS numbers: 31.25 .Eb, 36.10.Dr, 31.30.Jv, 3I.15.Pf, 02.70.-c 

Positronium ion (Ps~), consisting of two electrons and a positron, is the only known three-body bound state free from 
nucleons. Its existence was predicted by Wheeler in 1946 fT] and confirmed experimentally by Mills in 1981 f^. Only the 
ground state is stable against a dissociation into positronium and an electron (see [3] for an extensive review of its properties 
and references). Electron-positron annihilation limits the Ps^ lifetime to about half a nanosecond, as first reported in |4]. Here 
we determine relativistic and radiative corrections to the annihilation in a three-body bound state and predict the VsT decay rate 
with a 6 parts per million precision, 

r(Ps") = 2.087085(12)ns"^ (1) 

What makes the Ps^ ion particularly interesting is that its theory is very clean, albeit somewhat technically challenging. With 
a very good accuracy all but electromagnetic interactions in Ps^ can be neglected. Also the charge distribution of constituents 
is well known (point-like), unlike in atoms and ions containing nuclei. Quantum electrodynamics (QED) suffices to describe all 
properties of VsT . 

On the other hand, Ps^ is a three-body system and thus its wave function is not known analytically even in the non-relativistic 
approximation. This complicates theoretical investigations but also provides an opportunity to develop and test advanced com- 
putational techniques. Those new methods are important for other systems such as the hydrogen ion, the molecule H2 , and the 
helium atom. 
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FIG. 1: The main decay channel of Ps^ (a), and an example of a correction to it, the three-photon aimihilation (b). 

The e+e" annihilation proceeds fastest when the pair is in a spin-singlet state, like para-positronium (pPs), in which case two 
photons can be produced (see Fig. [It a))- If the pair is a spin triplet, like ortho-positronium (oPs), the decay results in an odd 
number of photons. Fig. Htb). Interestingly, unlike ortho-positronium, VsT can also decay into a single photon. However, this 
channel is very rare [5, 6]: all three constituents have to overlap to transfer momentum to the non-annihilating electron. The 
three-photon decay is much more likely, but still much slower than the spin-singlet two-photon process. 

The spatial wave function of VsT is symmetric with respect to the two electrons. For the total wave function to be antisym- 
metric, the two electrons must be in the spin-singlet state. It is convenient to think of Ps^ as consisting of a positronium core 
and a loosely bound electron |7]. This picture reveals the main features of the Ps^ lifetime. When e+ meets one of the e^, the 
odds are about one in four that their spins form a singlet. Thus Ps^ lives about four times longer than pPs. 

A variational determination of the VsT wave function |7] confirms this elegant argument. Furthermore, if this decay is so 
similar to that of pPs, the same O (a) corrections apply |8]. In the same order, also the three-photon annihilation must be 
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accounted for |0]. Together, this led to the theoretical prediction for the Ps decay width 10], 

rjf 3 ^ 2.086(6) ns-\ 

where the size of the O (a) corrections was used to estimate the uncertainty (ld\ (see also jU, 12, Tsll ). 

Recent measurement jlOtl agrees with this prediction and approaches its precision, Foxp = 2.089(15) ns^^. It is anticipated 
that the new intense source of positrons at the Garching reactor FRM-II will be used to decrease the experimental error by a 
factor of 4 or 5, below the uncertainty in Eq. (|2]i. 

Motivated by this effort, we undertook to improve the theoretical precision by determining all O (a^) effects. Ps^, a non- 
relativistic bound state, is well described by the Schrodinger equation. Its leading-order decay rate is 

ro-2^mea5((53(ri2)), (3) 

where me is the electron mass, ri2 is the distance between the positron and the electron which annihilates, and the mean value 
refers to the ground state Ps^ wave function 

i2T3 — T2i3 
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Throughout this paper we use l/amg, arrie, and a^nie as units of length, momentum, and energy (we also set c 
except in the last Eq. (l22l i). Thus, {5^{ri2)) in Eq. Q, as well as all mean values to follow, are dimensionless. 

Relativistic effects, spin of the electron, and short-distance exchanges of photons with a virtuality O (rrie) are not accounted 
for by the Schrodinger equation, which includes only Coulomb potentials among the three constituents. Like in other non- 
relativistic systems lll4ll . these additional effects can be treated as perturbations and organized in a series in a, 
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The first-order correction A, already discussed, includes corrections to the two- and the three-photon channels. 
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Some authors (e.g. 10]) hint at additional O (a) effects but in our opinion none other exist at this order. 

In the next order, four photons contribute fl^, corrections O (a) must be included in the three-photon decay 17], and 
O (a^) in the two-photon decay Iil8ul9t, i20i]. 



B = B-^'' + B^'< + B'^'^ . 
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The last term is the focus of this paper. It is a sum of several effects: square i?squaiod of the O {a) correction A^^; hard-photon 
corrections Shard to the e+e^ 77 process; and soft corrections to the annihilation amplitude Saa and the wave function B^^: 
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All corrections which affect only the annihilation amplitude have already been computed in the context of the pPs decay. Since 
they do not depend on the particular bound state, they apply to the present analysis without changes. 

The correction to the wave function i?wf , sensitive to the three-body dynamics, is the most challenging. As we will see below, 
it is divergent and cancels the divergence in Shard- The term — 2a^ In a in Eq. (|5]l is a remnant of those divergences. In order to 
regularize divergences, we work in = 3 — 2e spatial dimensions. Thus, the non-relativistic Coulomb Hamiltonian becomes 
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where pa, fab = fa — fi, are momenta and relative distances of the positron 1 and electrons 2, 3. 

The wave function correction arises due to relativistic effects, which are treated as a perturbation and described by the Breit 
Hamiltonian, 



^(4) ^ _y^p1^j2,^Ls^^,^^) + u 
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where the Pauli matrices are labeled with the number of the fermion they are acting on. It is convenient to evaluate separately 
the spin-independent part h['^\ and the spin-dependent part i?2^'', B„f = Bhi + Bh2 + ^■ 

The effect of this perturbation is the following replacement in the formula for the decay rate, Eq. (O, 



(<53(ri2))^2(j'^(ri2) 
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Here (^^^^^y is the Green's function of the lowest-order Schrodinger equation and the prime indicates the exclusion of the 
ground state. The appearance of divergences is the main obstacle in the evaluation of this correction. They originate from 
ri2 — > (ultraviolet limit), where the Breit Hamiltonian is not a valid description of the dynamics. Indeed, when one accounts 
for the hard photons, Eq. (|9]l, divergences cancel. 

In analogy with the earlier work on positronium and helium ll2lll22il23ll24ll . we rewrite the matrix element in Eq. ( fTSl l such 
that the divergences appear only in the coefficient of one operator, namely 6^{ri2). To this end, we rewrite ||24|1 the delta-function 
as 



4^5'^(ri2) = 47r<5'^(ri2) + \ Hq - E 
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This equation implicitly defines S'^, less singular than S''^. The most singular part is in the anticommutator in the second term. 

= / — l^*) (^'|, where / is the identity operator Hence, divergences 



This term cancels the Green's function, (E — Ho) „ 
appear only in first-order elements and are easier to extract. 
In the spin-independent part we find 



(17) 



from which we can determine the value of Bhi in terms of the ground-state mean values listed in Table U Among them, the 
regularized cubic operator is defined by Eq. (1.5) in Ref. 12411 . 

In the spin-dependent part, the effect of Pauli matrices in i?2a),, Eq. ( fT4l i. is evaluated with the spin wave function in Eq. (|4|l 
and represented by constants Aab for each pair of fermion lines: A12 = —2 — 6e, A13 ~ —A23 = —2. We keep e only in the 
coefficient of the divergent part. After this simplification of spins, we find 
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In terms of the operators in Table|I] using the symmetry r2 ^ fs ™d the virial identity 2E = {V), we get 



2tt ((5^(ri2)) Bh2 = t; i^^7 - "f^S - Vl4 - I'lS + ^'is) + 7 (^'9 + Wio) + ^^-^W20 + 7W2I + ^24 + 7^(«25 " 2u2o). (19) 
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The numerical values in TableUare obtained with a variational method. The trial wave function is expanded in a 1000-element 
set of exponential functions lEsll 
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TABLE I: Operators affecting the Ps wave function and their ground state mean values. We denote V12 = V + 1 /ri2. 



i Operator Oi Vi = (Oi) 



1 E'^/ri2 0.023 327 6 

2 ^^12/^12 0.033 945 

3 2Epi/rw -0.014 986 7 

4 -2EVi2/ri2 -0.015 844 8 

5 -2piV/ri2 0.090 907 3 

6 -P3/ri2 -0.029 850 5 

7 -47r53(ri3)/ri2 -0.051014 

8 47r5^(r23)/ri2 0.001 794 7 

9 3E/rl2 -0.219 554 9 

10 -3Vi2/rf2 0.022 092 3 

11 -(pi X P2)'(l/ri2)(pi X p2y -0.001 809 5 

12 -(4 7r/3)<5^(ri2)pi -0.003 350 2 

13 -pi {3ri2ri,/rl2 - 5'' /ri2) P^JI -0.000 275 1 

14 ri2 ■ ri3/(2 r?2r?3) -0.000 562 9 

15 ri2 • r23/(2 ri2rii) -0.001 413 3 

16 -Pi(l/ri2)P3 -0.037 118 4 

17 -pi(lAi2)P3 -0.008 951 8 

18 -EaPa(lA?2K/2 -0.077 005 7 

(5^Vr,, + ri, r^„,/rf ,) /ri2 pi 0.296 062 9 

20 27r5^(ri2) 0.130 270 5 

21 2P(l/r^) 0.014 113 8 

22 H^"' -0.072 738 1 

23 4 7r(53(ri2) 7— J— ^ J/f' 0.178 732 5 

24 Y.a<b 47ryl,ft3^(ri2) 7r(?3(r„,) 0.334 788 9 

25 27r5^(r23 ) 0.00 1 07 4 4 



where a, h, c are chosen randomly, with a homogeneous distribution, from two k sets defined by variational boundary con- 
ditions Aik < Ofci < A2fc, Bik < &fei < i?2fe, Cik < Cki < C2fc. Two (or more) sets allow one to match the be- 
havior of the wave function at various distance scales and improve accuracy. We found the non-relativistic energy value, 
E — —0.262 005 070 232 980(1), that agrees with an even more accurate earlier result ll26ll . Previously obtained mean values of 
^'^{rab) fi3,^2f], and non-singular products of l/vab ifisll are also confirmed. Finally, the mean value of the spin-independent 



part of the Breit Hamiltonian ijl^agrees with Ref . ||27|] . Crucial for the decay is the mean value of the delta-function, obtained 



using the representation of Ref. 112811 . 



(^3(ri2)) = 0.020 733198 0051(2). (21) 

This value agrees with the one found in II27I1 and somewhat improves its accuracy. 

For the new evaluation of the Ps~ decay rate we use a = 1/137.03599911(46) and the atomic unit of time a^nieC^/fi — 
10^^ s/2.418884326505(16) lH. Our final result in Eq. O is obtained using 

r(Ps-) = 27r^^(l + C) (SHn^)) , (22) 

where the correction C is given in Table HH and we use Eq. ( 1211 1. The last two corrections listed in Table lUrefer to the third order 
in a. The leading quadratic logarithm was found in 1, 3 0,1 and is valid for positronium atoms as well as for the ion. The linear log 
(the last correction) has not yet been calculated for Ps~ . However, it is known for pPs and oPs llsii [32i [33I1 . We expect its value 
for Ps^ to be close to that for pPs and use the latter as an estimate. We assign this correction a 100% uncertainty, which also 
conservatively estimates non-logarithmic higher-order effects 1,34.1 . 
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TABLE II: Corrections to the width of Ps 



Correction 


Value 




0.002 693 245 




-0.005 882 770 


-2a^ In a 


0.000 524 019 




0.000 001 480 




-0.000 064 352 


^ ^squared 


0.000 008 652 


„2 ofln 
" i'hard 


-0.000 218 3(34) 


a' Baa 


0.000 017 750 




0.000 078 366 




0.000 122 185 


3a^ln^ a/(27r) 


-0.000 004 491 


2.5(2.5)a=' In a 


-0.000 004 8(48) 


Total C 


-0.002 729 0(59) 



Note that Ref. ll27ll includes a prediction of r(Ps~ — > 77) with a seemingly higher precision than ours. That result, however, 
does not include any corrections beyond the tree level (this corresponds to setting C = in our Eq. ( |22] |) and its error estimate 
includes only the numerical uncertainty of the variational calculation in IztIi . 

Another experimentally interesting quantity is the branching ratio of the three-photon decay. We find 



BR(Ps" 777) = 



_ r(Ps-^777) 



r(Ps~ 
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A£''') - -A^'<a^ In ■ 



= 0.002 635 



(23) 



The uncertainty is due to the unknown O (a^) corrections to the decay Ps 111- Only the logarithmic term is known in this 
order [|35], and we take half of its value to estimate the uncertainty. 

The structure of corrections found in this study confirms the picture of Ps" as an electron loosely interacting with a positron- 
ium core 1 7 ] . The mean value in Eq. ( |2TI ) is very close to that obtained with a neutral positronium, neglecting the second electron, 
l/(167r) ~ 0.01989. Also, in the sum of all effects in Table HH there is a significant cancellation between the hard effects i?hard 
and the soft ones £?hi + i^H2 + i^aa, also observed in positronium ifioll . It would be interesting to understand the origin of this 
cancellation, which for now remains an open question. 

The accuracy we have obtained for the decay rate is 6 parts per milUon, about 500 times better than the previous best prediction, 
Eq. (O. Further progress in the theory of the Ps~ decay requires the logarithmic term O (a* In a) and improved hard corrections 
O (a^). However, the accuracy obtained in the present paper is sufficient for the foreseeable future. It exceeds the anticipated 
accuracy of Garching measurements by about a factor of 200. We are thus prepared for the new data and are looking forward to 
this intriguing test of three-body bound-state QED. 

We thank K. Pachucki and D. Habs for helpful discussions. This work was supported by the Alberta Energy Research Institute, 
under the COURSE Program, by Science and Engineering Research Canada, RFBR (grants # 06-02-16156 and # 06-02-04018) 
and DFG (grant GZ 436 RUS 1 13/769/0-2). 
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